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We formulate the expansion for the mass of the nucleon as a function of pion mass within chi- 
ral perturbation theory using a number of different ultra-violet regularisation schemes; including 
dimensional regularisation and various finite-ranged regulators. Leading and next-to-leading order 
non-analytic contributions are included through the standard one-loop Feynman graphs. In addition 
to the physical nucleon mass, the expansion is constrained by recent, extremely accurate, lattice 
QCD data obtained with two flavors of dynamical quarks. The extent to which different regulators 
can describe the chiral expansion is examined, while varying the range of quark mass over which 
the expansions are matched. Renormalised chiral expansion parameters are recovered from each 
regularisation prescription and compared. We find that the finite-range regulators produce con- 
sistent, model-independent results over a wide range of quark mass sufficient to solve the chiral 
extrapolation problem in lattice QCD. 



I. INTRODUCTION 

The importance of incorporating the chiral non- 
analytic behavior of hadronic observables in the extrap- 
olation of lattice QCD simulation results has been high- 
lighted in numerous studies over the past few years [1-9]. 
In an attempt to avoid model dependence, early lattice 
QCD extrapolations considered simple polynomial func- 
tions of the pion mass. However it is now widely accepted 
that the model-independent, chiral non- analytic behavior 
predicted by chiral perturbation theory (xPT) must be 
incorporated in any quark mass extrapolation function 
in order to preserve the properties of QCD. 

What is more controversial is the manner in which 
the effective field theory is regulated. Historically, most 
formulations of xPT are based on dimensional regulari- 
sation and surprisingly other regularisation schemes are 
often regarded as models. However, the physical predic- 
tions of the effective field theory must be regularisation 
and renormalisation scheme independent, such that other 
schemes are possible and may provide advantages over 
the traditional approach. Indeed, Donoghue et al. Q 
have already reported the improved convergence prop- 
erties of effective field theory formulated with what they 
call a "long-distance regulator," such as a dipole or Gaus- 
sian form, designed to suppress large three-momentum 
contributions to the loop integrals of heavy-baryon xPT. 

The origin of the controversy lies in the phenomeno- 
logical interpretation of the long-distance regulator as a 
form factor describing the finite size of the meson-cloud 
source. In this light, it is common to think of the reg- 
ulator as a model, removing the possibility of system- 
atic improvement through the calculation of higher order 
terms. However, the interpretation of the regulator as a 
physical form factor is valid only when the terms analytic 
in quark mass — the terms involving even powers of m w 
m Eq.(|) below - are either set to zero or determined 
independent of data by some model of hadron structure. 
While terms are not constrained by the effective field the- 
ory, one can maintain model-independence by including 



them as parameters fit to data, order by order in the ex- 
pansion. However there is an obvious corollary of this 
discussion. To suppress the size of the contributions of 
the analytic terms, one should consider the use of a phe- 
nomenologically motivated regulator. The suppression 
of the higher-order analytic terms in the corresponding 
residual expansion that may result from a good choice 
holds the promise of improved convergence properties in 
that expansion. 

It is extremely important to explore this possibility of 
improved convergence. At present it is clear that to link 
full dynamical-fermion lattice QCD simulation results for 
hadron masses with dimensionally- regulated xPT, new 
simulation results are needed in the range from the phys- 
ical pion mass to 0.1 GeV 2 |n[]; a daunting task even 
with the multi-Teraflops resources available over the next 
five years. Perhaps even more important is the formal 
issue of the formulation of 5?7(3)-flavor chiral perturba- 
tion theory, associated with a strange-quark mass, which 
likely lies beyond the applicable range of dimensionally- 
regulated chiral perturbation theory for most baryon ob- 
servables. 

This study will quantitatively compare the chiral ex- 
pansion of the nucleon mass for six different regulari- 
sation schemes, including traditional dimensional regu- 
larisation and a variety of long distance regulators. Our 
focus will be to determine over what range in quark mass 
(or squared pion mass) all chiral expansions agree at the 
1% level. It will become apparent that the long distance 
regulators do indeed have dramatically improved conver- 
gence properties and can imitate the dimensionally reg- 
ulated chiral expansion over a wider range than the con- 
verse. This range will also be reported. Finally, we will 
discuss the physics, explaining why dimensional regular- 
isation is a poor choice of regulator for effective field the- 
ory, in the context of lattice extrapolations. On the other 
hand, we reach the exciting conclusion that effective field 
theory formulated with a smooth long-distance regulator 
provides model-independent chiral expansions agreeing 
at the 1% level over a range exceeding < < 0.8 
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GeV 2 . Thus the chiral extrapolation problem has been 
solved. 



As mentioned above, coefficients of low-order non- 
analytic contributions are known, with ]15L |lq| 



II. CHIRAL EFFECTIVE FIELD THEORY 

Chiral perturbation theory (%PT) is a low-energy ef- 
fective field theory for QCD. Using this effective field 
theory, the low-energy properties of hadrons can be ex- 
panded about the limit of vanishing momenta and quark 
mass. In particular, in the context of the extrapolation of 
lattice data, xPT provides a functional form applicable as 
m q approaches zero or, equivalently, m n — > through the 
Gell-Mann-Oakes-Renner (GOR) relation m 2 cx m q |L2] . 
Goldstone boson loops play an especially important role 
in the theory as they give rise to non-analytic behaviour 
as a function of quark mass. The low-order, non-analytic 
contributions arise from the pole in the Goldstone boson 
propagator and hence are model-independent [fl3|| . How- 
ever, because the analytic variation of hadron properties 
is not constrained by chiral symmetry, the expansions 
contain free parameters which must be determined em- 
pirically by comparison with data. 

Effective field theory then tells us that the formal ex- 
pansion of the nucleon mass about the SU(2) chiral limit 
is: 

M N = a + a 2 m 2 + a^m\ 

+CT N7r (m 7r , A) + (7A7r(m-7r, A) + . . . , (1) 

where obtt is the self-energy arising from a Bit loop 
(with B = N or A) and A is a parameter associated 
with the regularisation. These N and A loops generate 
the leading and next-to-leading non-analytic (LNA and 
NLNA) behaviour, respectively. The expansion has de- 
liberately been written in this form to highlight that the 
theory is equivalently defined for an arbitrary regulator. 
This equivalence has been demonstrated by Donoghue 
et. al. [[Toll , where the convergence properties of SU(3) 
XPT were improved by implementing a finite-range reg- 
ulator (FRR), which they term "long-distance regulari- 
sation" . More recently this work has been extended to 
include decuplet baryons O]. 

The traditional approach within the literature is to 
use dimensional regularisation to evaluate the self-energy 
integrals. Under such a scheme the NNt: contribution 
simply becomes ON^m^, A) — > clnati^ and the analytic 
terms, a„m™ (with n even), undergo an infinite renor- 
malisation. The A contribution produces a logarithm, so 
that the complete series expansion of the nucleon mass 
about m T = is: 



M 



a = c + c 2 m 2 + c±m% 

+clna"4 + cnlna?t4 + 



(2) 



where the <Zj have been replaced by the renormalised (and 
finite) parameters C;. 
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Although strictly one should use values in the chiral 
limit, we take the experimental numbers with = 1.26, 
Jtc = 0.093 GeV, the nucleon-delta mass splitting, A = 
0.292 GeV and the mass scale associated with the loga- 
rithms will be taken to be 1 GeV (i.e. essentially 47r/ ff ). 

One would expect that it would be possible to reach 
the physical pion mass from the chiral limit with this ex- 
pansion truncated beyond c§m\. Unfortunately, the con- 
vergence of such an expansion seems to break down fairly 
quickly at higher pion masses. It is also important to re- 
alise that Eq. (^) was derived in the limit m^/A <C 1. At 
just twice the physical pion mass this ratio approaches 
unity. Mathematically the region m w rs A is dominated 
by a square root branch cut which starts at = A. Us- 
ing dimensional regularisation this takes the form [fL7| : 
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(4) 



for m-x < A, while for > A the first logarithm be- 
comes an arctangent. Clearly, to access the higher quark 
masses in the chiral expansion, currently of most rele- 
vance to lattice QCD, one requires a more sophisticated 
expression than that given by Eq. (0). 

Even ignoring the An cut for a short time, the formal 
expansion of the N — > N tt — > N self-energy integral, 
(TNtt, has been shown to have poor convergence proper- 
ties. Using a sharp, ultra-violet cut-off, Wright showed 
pH that the series expansion, truncated at 0(m^), di- 
verged for iTi-x > 0.4 GeV. This already indicated that the 
series expansion motivated by dimensional regularisation 
would have a slow rate of convergence. 

The main issue of the convergence of this truncated se- 
ries, Eq. (^|) , is linked to the formalism in which it is de- 
rived from the general form of Eq. (|l|) . The dimensionally 
regulated approach requires that the pion mass should re- 
main much lighter than every other mass scale involved 
in the problem. This requires that TOtt/A^b <C 1 and 
m w /A <C 1. An additional scale, mentioned briefly 
above, is set by the physical extent of the source of the 
pion field. This scale, which is of order R 
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responds to the transition between the rapid, non-linear 
variation required by chiral symmetry and the smooth, 
constituent- quark like mass behaviour observed in lat- 
tice simulations at larger quark mass |l9|]. An alterna- 
tive to dimensional regularisation would be to regulate 
Eq. (pi) with a finite ultra-violet cut-off (mass parame- 
ter A)J which physically corresponds to the fact that the 
source of the meson cloud is an extended structure. 
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In summary, while the low-energy effective field theory 
can be very useful in describing the quark mass depen- 
dence of hadron properties at very low values of m w , its 
utility in the context of lattice QCD appears to have been 
limited by the tendency in the literature to focus on a 
single type of regulator (i.e. dimensional regularisation) . 
We now investigate other possible regularisation schemes 
in order to see whether they are able to ameliorate the 
problem. 



III. CONSTRAINING THE CHIRAL 
EXPANSION 

Here we investigate the effective rate of convergence 
of the chiral expansions obtained using different func- 
tional forms for the regulator, together with an analy- 
sis of the dimensionally regulated approach. To analyse 
the merits of various regularisation schemes, at best, we 
would require exact knowledge of how the nucleon mass 
varies with quark mass. Having only one value for the 
experimental nucleon mass we cannot determine the pa- 
rameters that govern the quark mass dependence of M» 
without taking some information from alternate sources. 

Lattice QCD provides a non-perturbative method for 
studying the variation of Mn with m % . In principle, this 
allows one to fix the parameters of the chiral expansion 
using data obtained in simulations performed at vary- 
ing quark mass. Lattice simulations of full QCD are re- 
stricted to the use of relatively heavy quarks and hence 
it is not clear, a priori, whether the effective field theory 
expansion is capable of linking to even the lightest sim- 
ulated quark mass where ~ 500 MeV. An analysis of 
each regularisation scheme is necessary to determine the 
effective range of applicability. 

Taking as input the physical nucleon mass and the lat- 
est lattice QCD results of the CP-PACS Collaboration 
pot enables us to constrain an expression for Mn as a 
function of the quark mass. The lattice results have been 
obtained using improved gluon and quark actions on fine, 
large volume lattices with high statistics |2£|. In this work 
we concentrate on only those results with m sea = m V ai 
and the two largest values of f3 (i.e. the finest lattice 
spacings,[p9[ a ~ 0.09 - 0.13 fm). This ensures that the 
results obtained represent accurate estimates of the con- 
tinuum, infinite-volume theory at the simulated quark 
masses. The lattice data lies in the intermediate mass 
region, with m 2 between 0.3 and 0.7 GeV 2 . One may ask 
whether the effective field theory has any applicability at 
this mass scale. The following analysis will answer this 
question in a model-independent fashion. 

In order to remove the bias of choosing any particular 
regularisation scheme we allow each scheme to serve as 
a constraint curve for the other methods. In this way 
we generate six (one for each regularisation) different 
constraint curves that describe the quark mass depen- 
dence of Mn. The first case corresponds to the truncated 
power series obtained through the dimensionally regu- 



lated (DR) approach, Eq. (||) . We work to analytic order 
m^., which is necessary in order to counter the large, non- 
analytic contribution arising at order rn^Xogm^. The 
second procedure (labeled "BP") takes a similar form 
but the branch point obtained from the dimensionally 
regularised N — > An transition is retained in its full func- 
tional form. That is, the logarithm in Eq. ^ is replaced 
by the full expression, Eq. (|4|), which ensures the correct 
non-analytic structure where the logarithm converts to 
an arctangent above the branch point. We refer to this 
form as the dimensionally-regulated branch-point (BP) 
approach. Finally, we use four different functional forms 
for the finite-ranged, ultra-violet vertex regulator. These 
are namely the sharp-cutoff (SC), 6 (A — k); monopole 
(MON), A 2 /(A 2 + fc 2 ); dipole (DIP), A 4 /(A 2 + fc 2 ) 2 ; and 
Gaussian (GAU), exp(— fc 2 /A 2 ). It is commonly assumed 
in the literature that the functional form chosen for the 
finite-range regulator will introduce model dependence 
p5| and these rather different forms are chosen in order 
to test this hypothesis. 

Since, in the case of finite-range regularisation, we al- 
low the regulator parameter to be tuned to the lattice 
QCD data points, we do not require an analytic term at 
order m%. The use of a finite-ranged regulator will implic- 
itly include a term of order m\ , together with a string of 
higher order terms. Tuning the regulator parameter op- 
timises the efficiency of the one-loop chiral expansion. In 
particular, it naturally suppresses the large non-analytic 
contributions which would otherwise dominate at high 
pion masses. This is in contrast to the dimensional reg- 
ulated case, where the inclusion of the term in ra% is 
a practical necessity because of the large, non-analytic 



contribution at order 



In) 



Fixing the values of the non-analytic contributions to 
their model-independent values, we have four remaining 
parameters to be constrained for each of the regularisa- 
tion schemes. All of the regularisation schemes are con- 
strained simultaneously to fit the physical nucleon mass 
as well as the lattice QCD data. The resultant curves are 
displayed in Fig. [j]. To the naked eye all of the curves are 
very much in agreement with each other. All of them are 
able to give an accurate description of the lattice data 
and match the physical value of Mm- 

Now the question that we pose is: if any one of these 
curves were presumed to produce an exact description of 
the mass of the nucleon as a function of m n , how well 
could an alternate regularisation scheme match it? All 
forms have been based on the same, low-energy effective 
field theory and therefore will be equivalent in the limit 
m q — > 0. Hence, an additional question of more direct 
practical importance is: over what range can a partic- 
ular regularisation scheme match an alternate scheme? 
Within the radius of convergence, physical conclusions 
must be independent of regularisation and renormalisa- 
tion. 

We show the best fit parameters for our constraint 
curves, M{m„), in Table | It is worth recalling at this 
stage that the parameters listed in this table are bare 
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FIG. 1: Various regularisation schemes providing constraint 
curves for the variation of Mm with pion mass. The short 
dash curve corresponds to the simple dimensional regulari- 
sation scheme (DR) and the long dash curve to the more 
sophisticated dimensionally regulated approach (BP), which 
keeps the correct non- analytic structure at the Air branch 
point. The four finite-range regulators (solid curves) are in- 
distinguishable at this scale. Oscillations of the dimensionally 
regulated schemes (dashed curves) about the solid lines are 
apparent. 



Regulator 


ao 


ai 


a a 


ae 


A 


DR 


0.882 


3.82 


6.65 


-4.24 




BP 


0.825 


4.37 


9.72 


-2.77 




SC 


1.03 


1.12 


-0.292 




0.418 


MON 


1.56 


0.884 


-0.204 




0.496 


DIP 


1.20 


0.972 


-0.229 




0.785 


GAU 


1.12 


1.01 


-0.247 




0.616 



TABLE I: Bare fit parameters obtained for the constraint 
curves, M(m^). All quantities in appropriate powers of 
GeV. Regularisation schemes include dimensional regularisa- 
tion (DR) , dimensional regularisation maintaining the correct 
An branch-point (BP), sharp cut-off (SC), monopole (MON), 
dipole (DIP) and Gaussian (GAU). 



quantities and hence renormalisation scheme dependent. 
If one is to rigorously compare the parameters of the ef- 
fective field theory, the self-energy contributions need to 
be Taylor expanded about = in order to yield the 
renormalisation for each of the coefficients in the quark 
mass expansion about the chiral limit. We refer to the 
appendix for the details of this procedure. A compari- 
son of the resulting quark-mass expansion for each of the 
regularisation schemes is shown in Table 0. The most re- 
markable feature of Table || is the very close agreement 
between the values of the renormalised coefficients, espe- 
cially for the finite-range regularisation (FRR) schemes. 
For example, whereas the variation in ao between all four 
FRR schemes is 50%, the variation in Co is less than half 
of one percent. For a-i the corresponding figure is 30% 



Regulator Co C2 ca 

DR 0.882 3.82 6.65 

BP 0.885 3.64 8.50 

SC 0.894 3.09 13.5 

MON 0.898 2.80 23.6 

DIP 0.897 2.84 22.0 

GAU 0.897 2.87 20.7 



TABLE II: Renormalised chiral expansion parameters for the 
constraint curves, M(m n )- AH quantities in appropriate pow- 
ers of GeV. Regularisation schemes include dimensional reg- 
ularisation (DR), dimensional regularisation maintaining the 
correct Air branch-point (BP), sharp cut-off (SC), monopole 
(MON), dipole (DIP) and Gaussian (GAU). 



compared with 9% variation in ci- ff one excludes the 
less physical sharp cut-off (SC) regulator, the monopole, 
dipole and Gaussian results for C2 vary by only 2%. Fi- 
nally, for C4 the agreement is good for the latter three 
schemes. 

The comparison between a± and C4 is especially impor- 
tant in order to understand why the FRR schemes are so 
efficient. Whereas the renormalised coefficients are con- 
sistently very large for the three smooth FRR schemes, 
the bare coefficients of the residual expansion are a fac- 
tor of 100 smaller! That is, once one incorporates the 
effect of the finite size of the nucleon by smoothly sup- 
pressing pion loops for pion masses above 0.4-0.5 GeV, 
the residual series expansion has vastly improved conver- 
gence properties. 

In contrast, in order to fit the nucleon mass data over 
an extended range of pion mass, the dimensional regular- 
isation schemes require bare expansion coefficients which 
are much larger, 30 times larger in the case of 04. Still 
these coefficients are not large enough to reach the con- 
sistent values of the smooth FRR results reported in Ta- 
ble ||. The failure of this method is amply illustrated by 
the incorrect behaviour of M/v as soon as one applies it 
even slightly outside the fit region - c.f. the behaviour of 
DR and BP for m 2 . between 0.8 and 1.0 GeV 2 shown in 

Fig. g. 



IV. COMPARISON OF VARIOUS 
REGULARISATION SCHEMES 

With these given constraint curves, M (m^), we wish to 
test how well an alternative regularisation technique can 
reproduce the same curve. For any well-defined quantum 
field theory, all physical results should be independent of 
the regularisation and renormalisation schemes. By do- 
ing a one-to-one comparison over different ranges of pion 
mass we are able to determine the effective convergence 
range of each scheme. 

With one constraint curve, M(m 7r ), we have five alter- 



■5 



6000 
5000 
Q 4000 




0.4 

m w 2 (GeV 2 ) 

FIG. 2: RMS area between the dimensionally regulated con- 
straint curve and other regularisation schemes, x, plotted as 
a function of the curve fitting window (0, m^). 
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FIG. 5: RMS area between the monopole regulated constraint 
curve and other regularisation schemes, x, plotted as a func- 
tion of the curve fitting window (0, rxi^y). 
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FIG. 3: RMS area between the improved dimensionally reg- 
ulated constraint curve and other regularisation schemes, \, 
plotted as a function of the curve fitting window (0, m^). 



FIG. 6: RMS area between the dipole regulated constraint 
curve and other regularisation schemes, \, plotted as a func- 
tion of the curve fitting window (0, m w ). 
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FIG. 4: RMS area between the sharp cut-off regulated con- 
straint curve and other regularisation schemes, x, plotted as 
a function of the curve fitting window (0, m^y). 



FIG. 7: RMS area between the Gaussian regulated constraint 
curve and other regularisation schemes, \, plotted as a func- 
tion of the curve fitting window (0, rn^y). 
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nate curves, F(m w ), containing free parameters, which 
may be adjusted to fit this constraint. We fit each par- 
ticular regularisation scheme to our constraint curve over 
some window, € (0,tovi/). In doing so we define a 
measure, Xi which describes the root-mean-square area 
between the curves 



X 



dm% [F(m„) — M{m^)Y 



(5) 



In practical calculations the integral over is approx- 
imated by the Ricmann sum by dividing the region into 



N segments of size h 



/N, 



N 



dmlf(ml)^hJ2f(ih): 



(6) 




where h is fixed to 0.001 GeV 2 . For each upper limit, 
rriyy, the measure \ is minimised in the parameter space 
of the test function form, F(jn n ). 

All regularisation prescriptions have precisely the same 
structure in the limit — » 0. As our test window moves 
out to larger pion mass the variation of x wm describe 
the utility of each expansion. As an indicative scale, 



with mw = 0.5 GeV 2 a value of % = 700 MeV 2 means °= 



that the fit function is, on average, within 1 MeV of the 
constraint curve. Accurate reproduction of the expan- 
sion coefficients also serves to test the convergence over 
a given range. 

The variation of x as a function of the upper limit of 
the curve matching window, mw, is plotted in Figs. ^ 
to 0. The first feature that one notices is that, if ei- 
ther dimensionally-regulariscd scheme is used as the con- 
straint curve, all schemes are able to describe it very 
well for m 2 up to 0.4 GeV 2 . On the other hand, we 
see from Figs. ^ through |?] that neither DR nor BP are 
able to reproduce the behaviour of the FRR constraint 
curves beyond about 0.25 GeV 2 . However, all of the FRR 
schemes are able to reproduce the other FRR schemes 
over a much wider range - in fact, they reproduce them 
extremely well over the entire range of pion mass consid- 
ered (up to 0.8 GeV 2 ). Again, the coefficients in Table |l| 
help us to understand this observation: it is a simple 
consequence of the improved convergence properties of 
the residual series expansion once any reasonable FRR 
scheme is implemented. 



V. CHIRAL COEFFICIENTS 

We are concerned with the chiral expansion properties 
about the limit of vanishing quark mass. It is important 
to test how well the low-energy expansion parameters are 
reproduced as the curve-matching window is increased. 
Firstly, we take the dimcnsionally regularised curve as 
our constraint — i.e., we assume this is the most accu- 
rate description of the real world. We show in Fig. || 
the renormalised expansion parameters, Co, of the other 



m w 2 (GeV 2 ) 

FIG. 8: Recovered Co from various regularisation schemes con- 
strained to the dimensionally regularised (BP) curve. Co is 
plotted as a function of the curve-fitting window {Q,m\ v ). 
The dashed horizontal line indicates the maximum deviation 
of Co to within 1% of the physical nucleon mass. 
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m w 2 (GeV 2 ) 

FIG. 9: Recovered cq from various regularisation schemes con- 
strained to the dipole regularised curve. Co is plotted as a 
function of the curve-fitting window (Q,mw)- The dashed 
horizontal line indicates the maximum deviation of Co to 
within 1% of the physical nucleon mass. 



regulators constrained to the BP curve, as a function of 
the curve-fitting window. Figure ^ shows a similar plot 
for the case where the dipole regulator is taken as the 
constraint curve. 

It is clear from Fig. || that, whatever scheme is used to 
extract cq from the BP constraint curve, it is determined 
with an accuracy better than 1% as long as the fitting 
window is smaller than 0.7 GeV 2 . We define our 1% cut 
such that the contribution to the physical nucleon mass 
is less than 1% (e.g. in the case of C2, this requires that 
the error in the quantity C2m^/Mfq < 1% at the physical 
point). 

Conversely, Fig. ^ shows that the dimensional regular- 
isation schemes fail to yield the correct value of Co at the 
1% level once the fitting window to the dipole constraint 
curve extends beyond 0.5 GeV 2 . Yet all the FRR schemes 
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FIG. 10: Recovered C2 from various regularisation schemes 
constrained to the dimensionally regularised (BP) curve, ci 
is plotted as a function of the curve-fitting window (0, m^). 
The dashed horizontal line indicates the maximum deviation 
C2 can take such that the error in the quantity C2m^ is less 
than 1% of the physical nucleon mass at the physical pion 
mass. 
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FIG. 11: Recovered C2 from various regularisation schemes 
constrained to the dipole regularised curve. C2 is plotted as 
a function of the curve-fitting window CO. m%, ) . The dashed 
horizontal line is as described in Fig. hoi 



are accurate at a level better than half a percent, what- 
ever window is chosen — with the three smooth schemes 
accurate to 0.1%. 

A similar story applies to the coefficient C2, shown in 
All methods yield the correct value of 



Figs, hd and 11 



C2 based on the BP constraint curve within 1% out to 
a fit window of 0.7 GeV 2 . The breakdown at this scale 
should come as little surprise, since the highest lattice 
data point used to constrain the BP curve lies at around 
0.7 GeV 2 . 



When matching to the dipole constraint (Fig. 11), the 
dimensionally regularised cases break down much earlier, 
at around 0.4 GeV 2 . All three smooth FRR schemes 
reproduce this coefficient over the entire range to similar 
accuracy as found for cq. 



As we commented in Section IB , consistent results are 
found for the renormalised coefficients C4 over the entire 
fit range < m 2 < 0.8 GeV 2 , provided a smooth finite- 
range regulator is used to regulate the effective field the- 
ory. C4 is compromised in the dimensionally-regulated 
schemes as it attempts to simulate higher-order terms 
absent in the truncated expansion. 

It is clear from these exercises that one cannot hope to 
use a series expansion based upon dimensional regularisa- 
tion to analyse data over a window wider than 0.4 GeV 2 . 
That means one would need to have sufficient, accurate 
lattice QCD data in this mass region to fix four fitting 
parameters before one could hope to trust the chiral coef- 
ficients so obtained. Within this region, it would also be 
necessary to have data very near the physical pion mass. 
This is certainly beyond the likely computational capac- 
ity of current collaborations in the next 5-10 years. On 
the other hand, it is apparent that by using the improved 
convergence properties of the FRR schemes one can use 
lattice data in the region up to 0.8 GeV 2 or greater. This 
is a regime where we already have impressive data from 
CP-PACS fH and MILC f2f|. In particular, the FRR 
approach offers the ability to extract reliable fits where 
the low-mass region is excluded from the available data. 
The consideration of this practical application to the ex- 
trapolation of lattice data will be discussed further below. 

Indeed, these results suggest that given data from the 
next generation of (10 Teraflops) lattice QCD machines 
currently under construction, the FRR schemes should 
allow us to extract Co and C2, independent of any model, 
at the 1% level and C4 at the level of a few percent. 



VI. DISCUSSION AND CONCLUSIONS 

We have observed that all the regularisation schemes 
produce consistent results for cq and C2, with an accuracy 
sufficient to predict the nucleon mass at the 1% level, 
provided the expansion is limited to pion masses within 
the range < m 2 < 0.4 GeV 2 . The curves are very 
similar with an RMS area between pairs of curves less 
than 0.001 GeV 2 . 

Bowever, it is the dimensionally-regulated curves that 
are unable to accurately reproduce the results of the 
finite-range regulated curves beyond the range < m 2 < 
0.4 GeV 2 . As contributions from higher order terms not 
included in the truncated dimensionally-regulated expan- 
sion become important, the coefficients of the existing 
terms move away from their correct values in order to 
compensate for missing terms. 

On the contrary, the finite-range regulated curves are 
able to reproduce the predictions of dimensional regular- 
isation at the 1% level over a wider range approaching 
< m 2 < 0.7 GeV 2 . This enhanced range of the finite- 
range regulated curves is largely a result of the flexibility 
afforded by the presence of a separate regulation scale, 
A, which may be adjusted while preserving the excellent 
convergence properties of the finite-range regulated ef- 
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fective field theory. Hence, one can conclude that the 
applicable range of finite-range regulated effective field 
theory in which model independent conclusions may be 
drawn is < ml < 0.7 GeV 2 . This range is sufficient 
to solve the problem of formulating SU (3)-flavor chiral 
perturbation theory where the strange quark gives rise 
to a pseudoscalar kaon mass of 0.5 GeV. 

However, it is the predictions of dimensional regular- 
isation that prevent this range in pion mass from be- 
ing larger. Indeed the predictions of each of the smooth 
finite-range regulated formulations agree at an extraor- 
dinarily precise level over the entire range of pion mass 
considered, < m\ < 0.8 GeV 2 . This suggests that the 
more restricted regimes discussed above are associated 
with a fault in the choice of dimensional rcgularisation 
as a regulator for effective field theory in this context. 

It is easy to see why this may be. In dimensionally- 
rcgularised xPT, the mass of the pseudoscalar gov- 
erns the scale of physics contributing to the loop inte- 
grals of xPT. As the pseudoscalar mass increases, the 
dimensionally-regularised results become dominated by 
short distance physics; precisely the regime where effec- 
tive field theory fails because of the internal structure of 
the hadrons represented by the effective fields. Ignoring 
this physics by using point-like couplings leads to a badly 
divergent series expansion at higher order which, in turn, 
makes the method unsuitable for the lattice extrapolation 
problem in the foreseeable future. 

Fortunately this incorrect short-distance physics, in- 
cluded in dimensionally-regulated loop integrals, can be 
removed by a more appropriate choice of regulator and a 
corresponding shift of the parameters of the chiral La- 
grangian. In the region where these expansions both 
converge their predictions must agree as long as the se- 
ries are not truncated too soon. In a truncated ex- 
pansion, these incorrect contributions are not removed. 
In dimensionally-regularised xPT they become large for 
increasing m^, introducing a significant model depen- 
dence in the truncated chiral expansion and giving rise 
to the catastrophic failure already apparent in Fig. 1 
for m 2 > 0.8 GeV 2 . Convergence of the dimensionally- 
regulated expansion is slow as large errors, associated 
with short-distance physics in loop integrals (not sup- 
pressed in dimensionally-regulated xPT), must be re- 
moved by equally large analytic terms. 

The enhanced pion-mass range of the finite-range reg- 
ulated curves is a consequence of the effective resum- 
mation of the chiral expansion that arises via the reg- 
ulator. For example, a Taylor series expansion of the 
dipol e-reg ulated result for the N — ► Nir self energy of 
Eq. ( A13 ) in the Appendix reveals not only the stan- 



dard LNA contribution, but also new non-analytic 
contributions proportional to m\ and higher powers. As 
shown by Birse and McGovern, the coefficient of the m% 
term is related to the Goldberger-Treiman (GT) discrep- 
ancy [^7j j3(J . The coefficients of these new terms involve 
the dipole regulator parameter A which can be optimally 
constrained in the fitting procedure to reproduce the cor- 



rect non-analytic structure displayed in the data. (In- 
deed, the sign and magnitude of the coefficient of m\ 
that we find agrees with the value found by Birse and 
McGovern - although the latter has relatively large er- 
rors.) Such terms can only arise at higher order in the 
dimensionally-regulated expansion as there is no effective 
re-summation of the series. 

It is interesting to note that the appearance of higher- 
order non-analytic structure, such as the m% term dis- 
cussed above, does not occur for the special case of the 
sharp c ut-of f regulator. An expansion of the arctangent 
in Eq. (All) provides a constant followed by odd powers 
of in-n such that when multiplied by m\ only analytic 
terms appear. This unfortunate feature of the sharp- 
cutoff regulator explains why it lies between dimensional- 
rcgularisation and the smooth finite-range regulated re- 
sults. While analytic terms can be effectively re-summed, 
the sharp-cutoff regulator suffers the same restrictions 
as dimensional regularisation, in that higher-order non- 
analytic behavior is not effectively re-summed. 

The precise agreement between the smooth finite-range 
regulated results over the entire pion-mass range consid- 
ered, < to 2 < 0.8 GeV 2 , confirms that the shape of 
the regulator is irrelevant provided that the regulator pa- 
rameter is optimised by the fit to the data. An optimal 
regulator (perhaps guided by phenomenology) effectively 
re-sums the chiral expansion, encapsulating the essen- 
tial physics in the first few terms of the expansion. The 
approach is systematically improved by simply going to 
higher order in the chiral expansion, introducing addi- 
tional analytic terms as appropriate to maintain model 
independence. 

It is often argued that dimensionally regularisation 
works well in the meson sector. However, we know of 
no reason that the findings discussed here for the baryon 
sector should not apply equally well there. Improved con- 
vergence of finite-range regulated xPT for meson prop- 
erties will provide access to a wide range of quark masses 
and precision determinations of the low-energy constants. 

As discussed in the introduction it of great interest 
how to apply the features of effective field theory to the 
extrapolation of lattice QCD results. If one wishes to 
extract the low-energy constants from currently accessi- 
ble lattice simulations, it is necessary that any such pre- 
scription is not sensitive to the choice of regularisation or 
renormalisation. Here we describe how the FRR schemes 
achieve this result. 

With current lattice data typically restricted to m 2 > 
0.3 GeV 2 one must have a sufficiently wide window in 
order to accurately determine four fit parameters of the 
chiral expansion. Thus for the test of model-dependence 
in extrapolated results we presume perfect data up to our 
maximum considered range, m 2 = 0.8 GeV 2 . We fix the 
pseudo-data by the dipole regularised curve, the dipole 
being chosen as one of the three regulators which give 
best agreement over the widest range. Data points are 
chosen to be equally spaced by 0.05 GeV 2 , from 0.8 GeV 2 
down to some minimum available point at m\ . Fig- 
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FIG. 12: Extrapolated Co as function of lowest data point. 
Data points are spaced 0.05 GeV 2 apart and the maximum 
point is 0.8 GeV 2 . 
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Data points are spaced 0.05 GeV 2 apart and the maximum 
point is 0.8 GeV 2 . 
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APPENDIX A: RENORMALISATION 

Here we summarise the renormalisation prescription 
for various regularisation schemes. There are two types 
of diagrams which need to be considered in this analy- 
sis. The first of these is the NNw-type diagram, while 
the second corresponds to the case of a non-degenerate 
intermediate state, like the NAir. 

Within the heavy-baryon, non-relativistic approxima- 
tion the integral describing the N — > Nn diagram is writ- 
ten as 



dk 



\k) 



k 2 



(Al) 



where the factor u{k) is the Nit vertex regulator. The in- 
tegral has been defined in such a way that the coefficient 
of the m 3 contribution is normalised to unity. For the off- 
diagonal contribution the loop-integral can be expressed 



dk 



k A u 2 (k) 



\Jk 2 + to 2 (A + \fT? 



(A2) 



with A the A-N mass-splitting, defined as positive for a 
heavier intermediate state. The normalisation is chosen 
to agree with the previous equation in the limit A — > 0. 
The various functional forms for the finite-ranged regu- 
lators are listed in the text. 

Independent of the regulator, the Taylor expansion of 
the integral, /„■, behaves as 



ures [12| and |l3j show the reconstructed low-energy con- 
stants from the fitting window (m 2 ,0.8 GeV 2 ) with m 2 
in the range 0.4 > m\ > 0.05 GeV 2 . 

It is clear that the dimensionally regularised extrapo- 
lation schemes perform poorly. Even with 'perfect lattice 
QCD results down to 0.1 GeV 2 , cq is likely to be in er- 
ror by 7% and C2 by around 40%. On the other hand, 
one finds excellent agreement in the chiral parameters, Co 
and C2, obtained from the finite-ranged regulators, even 
with the lowest data point situated at 0.4 GeV 2 . This 
indicates that, using the technique of finite-range regu- 
larisation, one should be able to extract chiral coefficients 
at the 1% level from the lattice QCD data which will be 
obtained in the next few years. 



I n (m ~ 0) = &o + i>2m 2 + m 3 + 64771 4 + 



(A3) 



With the intermediate state being non-degenerate, the 
m 3 contribution becomes a logarithm and the expansion 
is given by 



I-K&{m ~ 0) = &oa + &2A"i 2 + 6 4 Af 
m 4 log m + . 



4 7T A 



(A4) 



With respect to the chiral expansion, Eq. ([!]), the self- 
energy contributions giving rise to the LNA and NLNA 
are given by 



o-N7r(nv,A) = C N /^(m^A) , 
^KjA) = C A /^(tott, A) , 



(A5) 



where for simplicity we have defined Cn = clna and 
Ca = — cnlna^kA/S. The renormalisation to obtain 
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TABLE III: Expansion coefficients for the finite regulators. 



the low-energy constants of the Taylor expansion, Eq. (||) , 
is then simply given by 



In Table HI we summarise the bi expansion coefficients. 
We were unable to obtain neat analytic expressions for 
the Gaussian regulator and hence all relevant calculations 
have been performed numerically. For both DR and BP 
schemes the bi of the 1^ diagram can trivially be treated 
as zero. The expansion terms, biA for the I,a diagram 
are more complicated and are summarised in the follow- 
ing. 

In the DR case, the 6^a can again be treated as zero. 
From the expansion of Eq. (|j), for the BP case one finds 
that the appropriate coefficients are given by (with renor- 
malisation scale 1.0 GeV) 



Co = a + Cn b + C\ &oa 
c 2 = a 2 + Cn b 2 + Ca &2A 
C4 = ci4 + Cn 64 + Ca &4A 



(A6) 



It is these renormalised coefficients, Cj, which have phys- 
ical significance and it is these that should be compared 
with phenomenological studies of chiral expansions. 



°0A 



h BP 
°2A 



h BP 
°4A 



log(4A 2 ) , 



2tt 
3 



16 7T A 



[l + 31og(4A 2 )] , 
[3 + 41og2 + 21og(A 2 )] . (A7) 
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For the sharp cut-off, expressions are given by 

1 



h sc 

°0A — 



°2A — 



3tt 



6A 2 A - 3AA + + 6A d log ■ 



1 



tt(A + A) 



A + A 

-A(3A + 2A) - 3A(A + A) log 



A 



-1 



hSC = 

4A 16^A(A + A) 2 



7 A 2 + 2 A A - 9 A 2 + 12 (A + A) 2 log 



A + A 
2 A A 



(A8) 



For the monopole we obtain 



.MON 
°0A 



lMON 
°2A 



A 4 



2tt(A 2 + A 2 ) 2 
A 2 

2tt(A 2 +A 2 ) 3 



A 2 

-2 A 3 - 2 A A 2 + 3 A 2 A tt + A 3 tt + 2 A 3 log ( 



A 5 + 6 A 3 A 2 + 5 A A 4 - 3 A 2 A 3 tt + A 5 ir 
A s 



-2 A A 2 (A 2 -3 A 2 ) log( =- 



lMON 
°4A 



13 A 8 + 54 A b A^ + 108 A 4 A 4 + 58 A z A e - 9 A w - 24 A 3 A 5 - 



16ttA(A 2 + A 2 ) 4 

+24 A A 7 tt - 6 A 4 (A 4 - 6 A 2 A 2 + A 4 ) log (4A 2 ) 

-6 (A 8 + 4 A 6 A 2 + 5 A 4 A 4 + 10 A 2 A 6 ) log (4A 2 ) 



(A9) 



For the dipole we obtain 



lDIP 
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A 4 



48tt(A 2 + A 2 ) 4 



I A 7 + 48 A 5 A 2 + 24 A 3 A 4 - 16 A A 6 - 3 A 6 A vr 



-27 A 4 A 3 tt + 27 A 2 A 5 tt + 3 A 7 vr + 96 A 3 A 4 log — 
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A 2 
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8 A y + 36 A 7 A 2 + 36 A A 4 + 188 A* A b + 180 A A 8 + 3 A b A 3 tt 



,5 A 4 



3 A6 



16 A 3 , 
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DIP 

a 4A 



+45A 4 A 5 tt- 135A 2 A 7 7r + 15A 9 7r-48 (5A 3 A 6 -3AA 8 ) log ^ 
-1 



54 A 12 + 326 A iU A 2 + 819 A" A 4 + 1062 A b A b + 1154 A 4 A 8 

48 A (A 2 + A 2 ) d tt 

+612 A 2 A 10 - 27 A 12 + 3 A 7 A 5 tt + 63 A 5 A 7 n - 315 A 3 A 9 tt + 105 A A 11 tt 
-6 A 8 (^35 A 4 - 42 A 2 A 2 + 3 A 4 ^) log (4A 2 ) - 6 A 2 ^3 A 10 + 18 A 8 A 2 



+45 A 6 A 4 + 60 A 4 A 6 + 10 A 2 A 8 + 60 A 10 log (4A 



(A10) 



af-i 2_/Y^ 2.A. 9 2 o 
1°^ = m z + -m 3 arctan 

3l IT TT 



Finally, we give the full expressions for the finite-range regulated integrals. For the simple diagram, /„., we obtain 

'A 

rMON = A 4 (2m + A) 
2(m + A) 2 ' 
DIP _ A 5 (m 2 + 4mA + A 2 ) 
ff ~ 16(m + A) 4 

For the more complicated integral, J,a, we obtain for the sharp cut-off 



%a - j- I 6A 2 A + 2A 3 - 6Am 2 - 3AAVA 2 + m 2 + 3A (-2A 2 + 3m 2 ) log( A + 

37T I ' m 
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For the monopole regulated integral we obtain 
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For the dipole regulated integral we obtain 
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ro = 0.5 fm p3[ 0]. This choice is ideal in the present 
context because the static quark potential is insensitive 
to chiral physics. 
[30] In dimensionally regulated x?T it arises only at the two- 
loop level, whereas in FRR it is already present at one 
loop due to resummation of the chiral expansion via the 



finite-ranged regulator. Furthermore, it is well known 
that a form factor at the nNN vertex with mass param- 
eters in the range shown in Table f| leads to a correction 
to the GT relation of the same size as the observed dis- 
crepancy. 



